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We consider two new aspects of Extensive Air Shower (EAS) development universality allowing to 
make accurate estimation of muon and electromagnetic (EM) shower contents in two independent 
ways. In the first case, to get muon (or EM) signal in water Cherenkov tanks or in scintillator 
detectors, it is enough to know the vertical depth of shower maximum X^ax and the total signal in 
the ground detector. In the second case, the EM signal can be calculated from the primary particle 
energy and the zenith angle. In both cases the parametrizations of muon and EM signals are almost 
independent on primary particle nature, energy and zenith angle. Implications of the considered 
properties for mass composition and hadronic interaction studies are briefly discussed. 

The present study is performed on 28000 of proton, oxygen and iron showers, generated with 
CORSIKA 6.735 for E'^ spectrum in the energy range log 10(£;)[eV] = 18.5 - 20 and uniformly 
distributed in cos^ 6 in zenith angle interval 6 = Q° — 65° for QGSJET II/Fluka interaction models. 

PACS numbers: 95.55.Vj, 95.85.Ry, 96.50.sd 



Introduction 



Mass composition of ultra-high-energy cosmic rays 
(UHECR) can be studied only indirectly with large EAS 
arrays. The contemporary measurement of longitudi- 
nal and lateral shower characteristics in hybrid experi- 
ments [ll-Q provides the possibility to combine several 
primary mass sensitive EAS parameters (such as depth 
of shower maximum and muon shower content) to achieve 
the best primary particle mass discrimination. Unfortu- 
nately, the lack of reliable information on hadronic in- 
teraction properties at these energies causes large un- 
certainties in the simulations of EAS characteristics and 
in turn brings large uncertainties in mass composition 
analysis results (see e.g. recent review [1]). For most 
of the interaction models the same experimental data 
show a lighter composition from Xmax and a heavier 
(generally with mass above iron) composition from muon 
shower content [sRsll . Recently proposed interaction 
model EPOS 0, llOj] seems to be the first model provid- 
ing quite consistent description of longitudinal and lat- 
eral EAS profiles due to increased muon production [ll[ . 
EPOS application to Hires-MIA data favors lightening of 
primary composition for energies above 10^* eV well in 
agreement with Hires recent results [H, [ll] (but com- 
pare with Yakutsk conclusions on mixed composition ob- 
tained with EPOS 1.60 for energies above lO^^ eV ^) 
and can partly help to reduce the 25% discrepancy be- 
tween Auger fluorescence and surface detectors energy 
scales [ll|. Nevertheless, this model has still to be thor- 
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oughly tested against various EAS data in the wide en- 
ergy range and, what's more, for energies E > 10^*, 
which are well beyond the accelerators reach, there is 
no way to judge if current EPOS (or any other model) 
properties correspond to actual nature of hadronic inter- 
actions. The difference between the various predictions 
of modern interaction models does not necessary cover 
all possible range of accelerator data extrapolations (see 
e.g. 15]) which is so wide that allows even to treat strong 
indications on primary mass increase for E > 10^^ eV, 
which can be derived from Auger data on ATmax and RMS 
of Xmax distribution [l2|, within pure primary proton 
fiux hypothesis by modification of the total cross sec- 
tion or by introduction of strong Feynman scaling viola- 
tion |16l - [l9| . Constraining of such interaction parameters 
freedom is crucial for resolving of UHECR mass composi- 
tion problem, but it seems to be hardly achievable on the 
basis of EAS studies alone and the data from the started 
LHC are considered of fundamental importance. 

In this p aper, u sing universality property of EAS de- 
velopment |20l - [25| . we propose two simple, independent 
and accurate methods to determine muon and EM shower 
contents and briefiy discuss a possible way to test and ad- 
just interaction models in a primary mass independent 
way. We also hope that the proposed EAS-universality- 
based correction of the interaction models will allow to 
perform mass composition analysis with the use muon 
EAS content in less interaction model dependent man- 
ner. 

The present study is performed making use of 28000 
showers, generated with CORSIKA 6.735 [265 for E'^ 
spectrum (and then weighted to E^^ spectrum) in the 
energy range log 10(i?) [eV] = 18.5 — 20 (with statistics of 
around 3000 showers for each primary in every of 3 energy 
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FIG. 1 EM and muon signals in showers from proton (red squares) and iron (blue crosses) in water Cherenkov tanks 
at 1000 m vs slant distance from shower maximum to the ground DG in log 10(i?)[eV] = 18.9 — 19.0 energy bin 
and 9 = 0°— 65° zenith angle range 
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FIG. 2 EM and muon signals in showers from proton (red squares) and iron (blue crosses) in water Cherenkov tanks 
at 1000 m vs vertical distance from shower maximum to the ground DG^ in log 10(i?)[eV] = 18.9 — 19.0 energy 
bin and 9 = 0° — 65° zenith angle range 



bins 18.5-19.0, 19.0-19.5 and 19.5-20.0) and uniformly 
distributed in cos^ 9 in zenith ang le interval 9 = 0° — 65° 
for QGSJET II d^-i^/Fluka [Sfl IsTj interaction models. 
The smaller set of around 10000 showers generated with 
EPOS 1.99/Fluka was also used for verification of the pre- 
sented universality properties. EM component thinning 
was set to 10^^, the observation level was at 870 g/cm^, 
geomagnetic field was set to Malargiie Auger site value. 
All longitudinal shower characteristics and charged parti- 
cles density D, which is effectively measured by detectors 
using scintillators [s^, were taken directly from COR- 
SIKA particle output files. The expected signal S in 
Cherenkov Auger-like tanks was calculated according to 
the sampling procedure described in [tI. IsSI with the use 
of the same GEANT 4 lookup tables as in Q • Differently 
from [7] in this work the muon signal Sf_i includes only 
signal from muons crossing the Cherenkov tank, while 
signal from EM particles, originating from muon decays. 



is included in the EM signal. 



I. SHOWERS AT THE SAME VERTICAL 
DEPTH OF MAXIMUM X^^^^ 

Of all aspects of universality of shower development we 
will be interested only in EM and muon signals depen- 
dence on the distance of shower maximum to the ground 
and on the zenith angle. Let's start from Auger-like ex- 
perimental setup and consider signal in water Cherenkov 
tanks at 1000 meters from the shower core. In this case 
the common way to express the universality of EM signal 
is to plot it against slant distance to the ground (Fig. [T] 
and also Fig. 1 in [1^), showing its quasi-independence 
from the nature of the primary particle. The muon signal 
functional dependence on slant distance to the ground 
DG is also very similar for both proton and iron, but 
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FIG. 3 Top; ratio of signals in water Cherenkov tanks Sp^/Sem at 1000 m vs vertical depth of shower maximum 
X^^^ in two energy bins. Protons — red squares, iron — blue crosses. Black line is the fit in the form ([T|) with 
parameters specified in the text. Bottom: distributions of relative difference between MC simulated muon signals in 
Cherenkov water tanks S*^^ and muon signals derived from the fit S"^* at 1000 m. Protons — red solid line, 
iron — blue dashed line. 



there is a shift in the normalization (Fig.[T]). Let us note, 
that since iron showers reach Xmax earlier than showers 
from protons, comparison of set of showers from p and Fe 
at equal DG means comparison between showers with dif- 
ferent zenith angles, but at the same development stage. 
Below, in this and in the next sections, we will consider 
other two possible geometrical situations: showers hav- 
ing the maximum at the same vertical depth ^^ax ^^'^ 
showers arriving at the same zenith angles. 

First, simply comparing shower characteristics depen- 
dence on vertical distance to ground DG^ , one finds a 
very interesting property (see Fig. [2|). In this case the 
similarity of the functional dependence of muon and EM 
signals on DG^ between p and Fe primaries is preserved, 
but now also EM signal normalizations are different. This 
happens because one compares showers, which have the 
same vertical distance from X^^^ to the ground, but pro- 
ton showers on average are more inclined than iron ones 
and therefore their EM component attenuates more on 
the way from X^nux to the ground. The ratio S'g^/S'gjy; 
turns out to be almost equal to the Sj^'^/S^ one and this 



allows to state a new shower universality property: the 
ratio of the muon signal to the EM one S'^/S'em is the 
same for all showers, reaching the maximum at the same 
vertical depth X^^^, independently on the primary par- 
ticle nature, primary energy and incident zenith angle (at 
least for the energy and angular ranges considered here). 
This property is illustrated in Fig. [3l where the depen- 
dence of Sf_i/SEM on Xj^ax P ^^"i primaries is shown 
in two different energy bins log 10(i?)[eV] = 18.5 — 18.6 
and 19.9 — 20.0. The functional dependence between 
X^ax ^^'^ Sfj^/SEM turns out to be very simple and quasi- 
universal for all energies and primaries. The function in 
the form 

X^^.^^AiSp/SEM + a)" (1) 

has been used to fit the data in 15 energy bins 
A log 10(i?)[eV]=0.1 and the fit parameters have been 
found to be stable across the entire energy range. Us- 
ing the functional dependence of X^^^ on (S^^/Sem) and 
Stot = Sem + one easily gets the equation, which al- 
lows to obtain the muon signal from shower vertical depth 
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FIG. 4 Means and RMS of distributions of relative 
difference between MC simulated muon signals in 
Cherenkov water tanks 5^*^ and muon signals derived 
from the fit Q S*^* at 1000 m (see also Fig.^, 
calculated with the unique set of parameters for all 
energy bins: A = 538, b = -0.25, a = -0.22. 
Protons — red, iron — blue. 



and total signal in water Cherenkov tanks: 



S\ 



■fit 



1 + l/((Xv^JA) 



i/b 



a) 



(2) 



We calculated the difference between the Monte-Carlo 
(MC) simulated muon signal and the muon signal 
obtained from the fit S'JJ* , examples of its distribution are 
shown in Fig. |31 In Fig 2] we plot behaviour of mean and 
RMS values of these distributions for various energy bins, 
obtained with the unique set of fit parameters A = 538, 
h = —0.25 and a = —0.22 which represent their aver- 
ages over 15 Alog 10(£')[eV]=0.1 energy bins. It is seen 
that the estimates of muon signals are unbiased with less 
than 1% deviation of the mean reconstructed muon signal 
from the MC one for all primaries and the RMS values are 
small: 8% for protons and around 5% for oxygen and iron 
(though we don't show results for oxygen which are al- 
ways between p and Fe, we used oxygen showers together 
with proton and iron ones to perform fits). Certainly, the 
application of specific coefhcients for each energy bin, or 
narrowing of zenith angle interval, or using of more so- 
phisticated fit functions can possibly slightly improve the 
performance of the method, but its simple and universal 
form already works very well. The described universal 
dependence of S^/S^m on X'^^^ holds true in the wide 
interval of distances (at the least from 200 to 1500 me- 
ters [S^l) from the core, though for distances closer to 
the core the function in the form ([Ij does not describe 
accurately the data in the entire angular range 0° — 65° 
and splitting the angular range in two regions or applying 
more complex parametrization may be needed. 

The same universality principle holds for the density 
of charged particles. We have performed the reconstruc- 



tion of muon densities using the dependence of the ra- 
tio (muon density D^)/ (electron density Df.) on X^^^ 
as in equation ([T]) with a unique set of parameters for 
all energy bins A = 475, b = -0.28, a = -0.09 (see 
Figs. [5l[6]). It is seen that when muons and electrons 
equally contribute to the detector signal, the shower fluc- 
tuations play a more important role and the accuracy of 
parametrization is only within 15%, though the estimate 
is still unbiased. 

In this context we would like to emphasize that since 
our study does not take into account any experimental 
error, we use the parametrizations only to demonstrate 
that the considered universality properties allow to per- 
form simple accurate estimates of the EAS characteris- 
tics and that the functional forms are really universal in 
terms of energy, primary particle nature and zenith angle. 
Certainly, their application to the real experimental con- 
ditions may require some corrections of the coefficients or 
even of the functional form of the parametrizations (es- 
pecially if core distances different from 1000 meters are 
considered). 

Hence, the S^/S-^m (Dfj/Dc) vs X^^^ universality al- 
lows to obtain accurate estimates of the muon signal 
(muon density) using simple parametrizations, which are 
almost independent on the primary particle nature, pri- 
mary energy and zenith angle for various types of ground 
detectors. Taking into account that the shower univer- 
sality pr operty was established for different interaction 
models [20l425l | . the proposed approach to muon con- 
tent derivation should not be specific only to QGSJET II 
and our test calculations with EPOS 1.99 confirm it (see 
Fig. [21). In case of EPOS 1.99 the fit in the form ^ 
still provides good description of the simulated data, but, 
as expected, fit coefficient are different from those for 
QGSJET II. 



II. SHOWERS AT THE SAME ZENITH ANGLES 

Another universality property follows from the study 
on showers arriving at the same zenith angles. In this 
case the average iron shower has to cross larger slant 
distance from ATmax to the ground with respect to the 
average proton shower and this almost equalizes EM 
signals (densities) for both primaries at the observa- 
tion level in the wide range of zenith angles (Fig. 
For the signal at 1000 meters in the Cherenkov water 
tanks, notable discrepancies between p and Fe showers 
EM components are observed for nearly vertical show- 
ers {6 < 18°, cos^{9) > 0.9) and very inclined ones 
{9 > 63°, cos^{e) < 0.2). In the first case the situa- 
tion is similar to the one considered in Fig. [TJ the path 
from Ainax to the ground for p and Fe showers is almost 
the same. For inclined showers the difference is caused 
by the EM halo from muon decays and larger number of 
muons in iron showers brings to a larger EM halo signal. 
We normalize S'em by primary energy to cancel the al- 
most linear growth of the signal with the energy in order 
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FIG. 5 Top; ratio of muon density to the electron one at 1000 m vs vertical depth of shower maximum ^^ax two 
energy bins. Protons — red squares, iron — blue crosses. Black line is the fit in the form ^ with parameters 
specified in the text. Bottom: distributions of relative difference between MC simulated muon density Z?^^ and 
muon density derived from the fit D^^ at 1000 m. Protons — red solid line, iron — blue dashed line. 



to obtain a universal description of S'em for all energies. 

Looking at the showers at different zenith angles one 
samples longitudinal showers profiles, for this reason it 
is natural to try to describe the dependence of the EM 
signal on cos(6') with Gaisser-Hillas type function, using 
cos(0) as variable instead of Xmax: 



S'em (cos (6'))/£; [VEM/EeV] = S\ 



■0 

EM 



X exp 



Cl 



cos(6') — Co 

Cl - Co 

cos(6l) 



A 



(3) 



where a 



(ci - co)/A; 



Co and ci (cosine of angle at which 5em='S'e]vi 



•^EM (signal at maximum), 
are fit 

parameters. The fit parameters S^yi and ci change 
by less than 10% and 3% correspondingly across en- 
tire range of energies (when one makes fits in 15 energy 
bins AloglO(£;)[eV]=0.1 from log 10(i;)[eV] = 18.5 to 
log 10(£')[eV] = 20.0), while co changes quite chaotically 
from to —20 (this causes to change also A). We have 
found that, fixing co (similarly to [7|) to any negative 



value within this range, we obtain a good universal fit 
and A changes in this case by less than 15%. Finally, 
we used the following average values (except co that was 
fixed to —3) of the coefficients S'^j^ — 2.53, co = —3, 
Cl = 0.96, A = 0.012. The resuhs of the fit and the dif- 
ference between the MC simulated EM signal and 
the EM signal obtained from the fit S^\^ are shown in 
Fig. |9l The accuracy of the reconstruction for all en- 
ergy bins is shown in Fig. [TU] and it is seen that one 
gets an unbiased estimate of S'em with RMS below 15% 
for proton and 13% for iron showers. Let us note that 
RMS of the distributions of relative difference between 
fit and MC EM signals reflects shower-to-shower fluctu- 
ations and that the application of fits with coefficients 
calculated accurately for each energy bin does not pro- 
duce any reduction of the RMS. 

The same fit procedure was also applied to the density 
of electrons D^- In this case good agreement between 
p and Fe EAS electron densities was found in the wider 
angular range — 63° and cq was fixed to zero, that gave 
the following set of fit parameters: 13° = 0.43 (instead of 
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FIG. 6 Means and RMS of distributions of relative 
difference between MC simulated muon density Djf ^ 
and muon density derived from the fit D^^ at 1000 m 
(see also Fig. [5]), calculated with the unique set of 
parameters for all energy bins A — 475, b = —0.28, 
a — —0.09. Protons — red, iron — blue. 



S'em)' Co = 0, ci = 0.96, A 0.069. The results of the 
fit and electron density reconstruction accuracy are pre- 
sented in Figs. [lll[T2l It is seen, that shower-to-shower 
fluctuations here are larger compared to the EM signal in 
water Cherenkov tanks and while the density estimates 
are still unbiased, their spread reaches 25% and 22% on 
average for proton and iron showers correspondingly. 

Finally, our test calculations presented in Fig. [T3l 
demonstrate that the universality of EM signal (density) 
dependence on zenith angle holds true also in case of 
EPOS 1.99. 



Conclusions 



forward, but any their modifications can be easily ac- 
complished and we believe that these universality prop- 
erties can be used in hybrid experiments for mass com- 
position studies, primary and missing energy estimates 
and for hadronic interaction model tests. In particular, 
we would like to dwell here on the problem of muon ex- 
cess in the real data compared to predictions of most of 
the interaction models 5-7]. Since the muon content of 
EAS is highly mo del- dependent (see e.g. [Ill [HI) and the 
UHECR mass composition is still unknown, this muon 
excess can be expressed only in terms of a relative excess 
with respect to the prediction of a given hadronic inter- 
action model for a given primary. In our point of view 
primary mass-insensitivity of the presented universality 
properties can be used to reveal inconsistencies in the 
hadronic interaction models without bias coming from 
the unknown mass of the UHECR. One of the possible 
strategies lies in simultaneous application of both consid- 
ered universality properties to the data in order to reveal 
hadronic model inconsistencies and, after taking them 
into account, to get concording estimates of muon and 
EM contents. Once it will be possible to get S'^x/S'em vs 
^max dependence for the real data (and if the universal- 
ity ([ij will be observed also there) it will be clear how one 
should rescale number of muons predicted by QGSJET II 
or by another model and this correction will be indepen- 
dent on unknown primary mass. Another possible appli- 
cation of the property ([T]) for muon content derivation in 
almost interaction model independent way with its con- 
sequent use for primary mass composition analysis will 
be presented elsewhere [s^. And finally, for hybrid de- 
tectors equipped with muon counters probably it will be 
possible to use 5'^/5em vs X^^^ universality for deter- 
mination of depth of shower maximum taking advantage 
of 100% ground array duty cycle in respect to 10% one 
of the fluorescence telescopes, but certainly this problem 
requires dedicated study. 



We have proposed two new EAS universality properties 
providing two independent ways to access EM and muon 
shower contents. We have shown that these properties 
can be described with simple parametrizations which are 
vahd in the wide energy (log 10{E)[eV] = 18.5-20.0) and 
zenith angle ranges, and independent on the primary par- 
ticle nature. Certainly, application of these parametriza- 
tions to the real experimental conditions is not straight- 
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water tanks -S*^ and muon signals derived from the fit at 1000 m. Protons — red solid line, iron — blue 
dashed line. 
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FIG. 10 Means and RMS of distributions of relative 
difference between MC simulated EM signals in 
Cherenkov water tanks and EM signals derived 
from the fit (O S'I'm at 1000 m (see also Fig.^, 
calculated with the unique set of parameters for all 
energy bins: 5}^^ = 2-53, cq = -3, ci = 0.96, A = 0.012. 
Protons — red, iron — blue. 
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FIG. 11 Top: electron densities at 1000 m vs cos(6') in two energy bins. Black line is the fit in the form (jSj with 
parameters specified in the text. Protons — red squares, iron — blue crosses. Bottom: distributions of relative 
difference between MC simulated electron densities D^'^ and electron densities derived from the fit at 1000 m. 
Protons — red solid line, iron — blue dashed line. 
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FIG. 12 Means and RMS of distributions of relative 
difference between MC simulated electron densities 
D^*^ and electron densities derived from the fit ^ D^*" 
at 1000 m (see also Fig. [11]), calculated with the unique 
set of parameters for all energy bins: = 0.43, cq = 0, 
ci — 0.96, A = 0.069. Protons — red, iron — blue. 




FIG. 13 EM signals in water Cherenkov tanks (left) and electron densities (right) at 1000 m vs cos(6') for EPOS 1.99 
in log 10(i?)[eV] = 18.8 — 19.0 energy bin. Lines show the fits in the form (|3]), but coefficients certainly differ from 
those for QGSJET II. The other designations are the same as in Fig. [SI 



